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USE OF MEAN ELEMENTS TO CALCULATE THE 
POSITION OF THE SUN, MOON AND EARTH 


by 


R. J. Sandifer 



the Earth, Moon and Sun are presented in tabular form for the 
epoch 12 hours ephemeris time, January 0, 1900. The elements 
are listed in six combinations of units: Three angular units of 
degrees, radians and revolutions and two time units of days and 
Julian centuries. 


A method is described for updating the numerical coefficients 
to any arbitrary epoch in order to eliminate double precision 
arithmetic in computer usage of these expressions. An algorithm 
is presented for the calculation of the selenocentric position vectors 
of the Earth and Sun neglecting periodic perturbations. 
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NOTATION 
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O) 
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C 
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u 

r 
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Superscript — 


none = 
Subscript — 


o — 
3 = 

none = 


Orbit Parameters 

semi-major axis 
eccentricity 

inclination of orbit plane 
longitude of ascending node 
argument of pericenter 
mean anomaly 
eccentric anomaly 
true anomaly 
equation of the center 
mean argument of latitude 
true argument of latitude 
magnitude of radius vector 

mean loi^tude, ecliptic and mean equinox of date 

mean longitude of pericenter, ecliptic and mean equinox of date 

Subscripts and Superscripts 
(When used with notation for orbit parameters) 

indicates center of coordinate system 

geocentric 

selenocentric 

heliocentric 

indicates body whose motion is being described by the orbit 
parameters 

Sun 

Earth 

Moon 


V 


! 


r 

8 ' 

O 

8 ' 


NOTATION 


mean obliquity of date 

inclination of mean lunar equator to the ecliptic 
Greenwich mean sidereal time 
Lunar mean sideral time 


1 


0 


0 


R,(a) 


0 cos a sin a 

0 -sin a cos a 


R2(/^) 


COS /3 0 -sin /S 

0 1 0 

sin /3 0 cos /? 




cos y sin 7 

-sin 7 cos 7 

0 0 


0 

0 
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USE OF MEAN ELEMENTS TO CALCULATE THE POSITION 
OF THE SUN, MOON AND EARTH 


INTRODUCTION 

The usage of the position vectors of the Sun, Moon, and Earth in digital com- 
puter programs designed to solve problems related to artificial satellites may 
be roughly divided into three categories: 

1. Orbital 

a. Position of central body 

b. Position of disturbing bodies 

2. Aspect and Orientation 

3. Shadow and Occultation 

The accuracy requirements placed upon the computer program in the calculation 
of these position vectors will usually depend upon two factors: 

1. Accuracy of Observations 

2. Experimental Requirements 

The method most commonly employed in the calculation of the position vectors 
of these bodies consists of a two step process. First, it is necessary to prepare 
a magnetic tape containing tables of the positions of the bodies as a function of 
the independent variable, time. The positions are usually tabulated at fixed in- 
tervals of time, the time interval being dependent upon the rates of change of the 
position of the particular body. The position of the Moon may require position 
entries every hovir, the position of the Sun every 24 hours, the position of the 
outer planets every 10 days, etc. The magnetic tape containing such information 
is referred to as an ephemeris tape. Secondly, to utilize this ephemeris tape, 
it is necessary to prepare two separate sets of computer instructions. The first 
set of instructions physically positions the magnetic tape to an appropriate initial 
entry and then reads off a number of tabulated entries from the tables. The 
second set of instructions interpolates between the tabulated entries to obtain 
the positions of the bodies for the desired calculation time. 
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The ephemeris tape so prepared usually consists of values of the position of 
the body to the highest available accuracy. In most cases, this accuracy is much 
greater than that required by experimental or observational considerations. 

This method has only one feature in its favor — it does not require any thinking 
on the part of the user. The user usually justifies the use of the method by estab- 
lishing a requirement for the most accurate available positions. In most cases 
this requirement will not stand up under a careful analysis of the problem. The 
disadvantages of the method include the usage of every known bad computer 
practice, not the least of which is the use of an input/output device to supply 
data which can be calculated by the computer. 

An alternate method is to design subprograms for the computer which will 
supply the required position vectors by direct calculation from available solutions 
to the equations of motion of the desired body.^ These solutions usually have the 
attractive feature of simplicity in two respects: 

1. Only one variable, time, is required to solve for the necessary elements 
of the orbit in order to obtain a position vector of the body. 

2. The solution for any one of the orbital elements is usually of the form 

e = + Cj t + + periodic terms. 

The following section will investigate various means of implementing this 
alternate method with emphasis upon the use of the "mean elements" of the orbit 
to obtain the position of the body. 


ALGORITHM FOR TWO BODY ELLIPTIC MOTION 

From the point of view of computer usage, the most useful solutions to the 
equations of motion of the Earth, Moon, and planets start on the assumption of 
"two body motion." 

To start the solution, it is assumed that there are no forces acting on the 
two bodies (e.g. the Earth in orbit around the Sun) which cause a deviation from 
"two body motion". On the assumption of the existence of the two bodies only, 
either as point masses or as homogeneous spherical masses, plus know initial 

Whis procedure was suggested in 1961 by R.L. Duncombe, currently Director of the Nautical 
Almanac Office, U.S. Naval Observatory (see page 453 “Transactions of the International 
Astronomical Union, Vol Xl-B, Proceedings of the Eleventh General Assembly, Berkeley 
1961")- To the best of the author’s knowledge, this approach has not yet gained widespread 
acceptance despite the many advantages it offers over other methods. 
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conditions in the form of a position and velocity vector of one body relative to 
the other at a specific reference time, t^j, or their equivalent, it can be proven in 
a strai^tforward manner that the orbit of one body relative to the second body 
describes a conic section. The position and velocity vector at any other time 
can be computed rigorously in terms of a set of six constant elements corresponding 
to the reference time and the time elapsed from . These six constant ele- 
ments can be computed rigorously from the initial conditions. In the event that 
the relative conic section is an ellipse, the corresponding constant elements are 
called "elliptic elements". One useful set of elliptic elements is the following 
set: 

a = Semi-major axis 
e = Eccentricity 

= Mean anomaly at time 
i = Inclination 
n = Longitude of node 
CO = Argument of pericenter. 

The elements a and e describe the size and shape of the ellipse and go de- 
scribes the position of the second bo(fy at the epoch time relative to the first 
or central body which is located at one focus of the ellipse. The three elements 
i , Q , o) describe the orientation of the ellipse relative to an arbitrary rectangular 
coordinate system centered at the central bo(fy. For two body motion all of the 
above elements are constant in time. These elements, or any equivalent set of 
six independent elements, can be used to find the position and velocity of the 
second body relative to the first body by means of the following algorithm for 
elliptic motion. 


Algorithm for Elliptic Motion 


1. Compute the mean angular motion 


n = 



where 

GM 

G = universal gravitational constant 
M = sum of the masses of the two bodies. 
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2. Compute the mean anomaly for any arbitrary time t 


g = go + n (t - tg) 

3. Compute the auxiliary angle E (eccentric anomaly) iteratively by means 
of the transcendental equation (Kepler’s equation) 

E - e sin E = g 

4. Compute position and velocity components relative to a rectangular co- 
ordinate system in the orbit plane, centered at the central body, with 
first principal axis directed towards pericenter, the second principal 
axis in the orbit plane, 90° away from the first principal axis in the 
direction of the motion of the second body, and the third principal axis 
normal to the orbit plane and directed in the right-handed sense: 





"a (cos E - e) 

q = 

^2 


a /l - sin E 




0 


r = a (1 - e cos E) 


/^i\ 




i//7a~ 


j - sin E ^ 
/l - cos E 
0 I 


where r is the magnitude of the radius vector. 

5. Rotate this position and velocity vector from the coordinate system de- 
scribed in 4 to the arbitrary cartesian coordinate system. 


r =R 3 (-n) (-i) R 3 (-OJ) q 

= Rq 
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where 


R = R3(-fi)Rj(-i)R3 (-a;) 


and 


r = Rq 

MEAN ELEMENTS 

If perturbations due to third bodies, non-sphericity of the two bodies, aero- 
d3Tiamic drag, or other causes are now taken into consideration, the above assump- 
tion of elliptic elements of constant value are no longer valid. If it is desired 
to have a solution of the equations of motion in terms of "elliptic elements" 
which will yield the position and velocity of the second body relative to the 
first body by means of the above formulae for elliptic motion, the desired set 
of "elliptic elements" are no longer constant in time. Letting the symbol e de- 
note in general a , e , g^j , i , fl , o) or any other equivalent "elliptic element," 
the solution now takes the form: 


e = d + Cj d^ + 63 d^ + P 

where P denotes periodic terms, ®o ’ constant numerical 

coefficients, and d is the time elapsed from , the epoch or reference time. 

This can also be written 


e = e + P 

m 

where d + d^ + 63 d^. is commonly called a mean element^ 

and e is called an osculating element, or an element of the osculating ellipse. 

The geometrical significance of osculating elements and the osculating 
ellipse may be described as follows. Assume that at a given time t all forces 
acting upon the two bodies in such a manner as to cause deviation from "two-body 
motion" are removed. From that time on, the two bodies will then describe 
simple elliptic motion. The elements which would then define this elliptic motion 
are called osculating elements. However, since perturbations are continually 


^There are various other definitions of the term “mean element* according to the method of 
solution of the equations of motion. The one described here is the most straightforward 
and useful one for computer applications. 
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acting on the two bodies, this ellipse, or osculating ellipse, is continually changing 
shape, size, and orientation in space. From a computational point of view, the 
important point is that if we have some solution which will yield a set of osculating 
elements, e = e (t) , for any time t, the position and velocity of the second body 
relative to the first may be easily found for that instant of time by means of the 
above algorithm for elliptic motion. 

Referring again to the general e3q>ression for the calculation of osculating 
elements 


e = + 6j d + + P 


there are five items worthy of note: 

1. If there are no perturbations, i.e. if we have simple two body motion, 
then 

e = or 6j = = 0 and P = 0. 

2. The expression for P, the periodic terms, usually takes the form of a sine 
and cosine series with small numerical coefficients. The arguments of 
the sine and cosine terms in turn usually consist of mean elements, also, 
so that the entire expression for e may be evaluated with only one vari- 
able, d, known, i.e. e = e(d) only. 

3. Since the coefficients of the individual terms of P are usually very small, 
the osculating value of any orbital element may be approximated to a high 
degree of accuracy by its mean element.^ 

e = = +e^d + e^d^ 


^One exception to this is the periodic term called evection in the calculation of the longitude 
of the Moon. The evection has an amplitude of approximately 1?3 and a period of about 32 
days. Acknowledgments are due to Mr. L. Werner of IBM who brought this to the attention 
of the author as a result of hi s computer program and discussions with Dr. P. Musen. 

^On page 113 of the “Explanatory Supplement to the Ephemeris”, an estimate of an accuracy 1 ' in 
the position of the inner planets is given using their mean elements. 

^Clemence^ G.M., “On the Elements of Jupiter”, Astronomical Journal November, 1946, Vol. 52 
Number 4, page 89, 
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4. With proper analysis of the problem, prerequisite accuracies may be 
obtained by inspection of the coefficients of the periodic terms of the 
elements, retainii^ only a sufficient number of the larger coefficients 
to meet the requirement of the problem. 

5. It is entirely feasible with current digital computer configurations and 
storage capacities to store the entire solutions of the equations of motion 
for the Earth, Moon, and planets in the memory of the computer. If 
these solutions are programmed in the form of a subroutine, the usage 
becomes extremely simple since the only input requirement is time, and 
the only limitation on the value of the time used is the limitation imposed 
by the solution itself. Any storage or computing time problems which 
may occasionally arise can be overcome by very simple programming 
techniques.^ 


^An interesting compromise between the practice of having ephemeris on tope and computing 
internally from the equations of motion has been developed for the planets by Mr. Lloyd 
Carpenter of the Theoretical Division. Mr. Carpenter’s technique consists of fitting a 
Chebyshev polynominal over discrete time periods. 


7 



METHOD OF UPDATING THE EPOCH OF MEAN ORBITAL 
ELEMENTS IN ORDER TO AVOID THE USE OF DOUBLE 
PRECISION ARITHMETIC 


The general expression for the calculation of the mean orbital elements of 
the Earth, Moon, and planets has the form 


e - e + e,d + e,d^ + e,d^ 

m o 1 2 3 


where is the numerical value of the mean element on any given Julian Date, JD; e^, 
Gj, Cj, and are numerical coefficients valid only for a specific epoch Julian 
Date, JD^; and d = JD - JD^. In the enclosed tables, the coefficients are valid 
for JD^ = 2415020. 0, the Julian Date equivalent to the calendar date, 12 hours 
ephemeris time, January 0, 1900. In order to use these formulae for the cal- 
culation of mean orbital elements during the 1960's, double precision arithmetic 
would be required in the calculation of the e ^d term in order to preserve ac- 
curacy. Extensive use of double precision arithmetic may be avoided by updat- 
ing the coefficients to be valid for a new epoch Julian Date, JD^, equivalent to 
some convenient calendar date. The term "calendar date" is used in the sense 
of year, month, day, and time of day. 

The required transformation is performed in the following manner. 


h 


JD, 


JD 


JDo 

old epoch (1900) 


1 

JD' 

O 

new epoch 

(arbitrary-to be chosen 
for convenience and ac- 
curacy of calculations) 


JD 

date for which 
mean element 
is to be calcu- 
lated 


2415020. 0 = Julian Date corresponding to calendar date, 12 hours 
ephemeris time, January 0, 1900 
Julian Date of new epoch 

Julian Date for which element is to be calculated 
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= JD^ “ = number of days from old epoch to new epoch 

d' = JD - JD^ = " " " " new " *' calculation date 

d = JD - JD = d + d' = " " " " old " " " ” 

o o 

An equivalent expression 


e - e' + el + e'd'^ + eld’^ 

m o ‘ 2 3 

where d ' is the numbpr of days elapsed from the new epoch JD^ , can be derived 
from the original expression by making the substitution d = d^ + d ' . 

+ e,(d^ + d') + 02(d^ + d')^ + 03 (d^ + d')^ 

= + e,d„ + e^d- + e3(d^" + d’^ + 2d„d') + e3(df + d'^ + 3dfd' + 3d„d’^) 

" (®o + + ^2^10 + ^3^l) ^ (^1 ^ ^ejd^ + 3e3df)d' + (ej + 3e3djd'^ + e3d'^ 

e; + e;d' + e-d*' + e'd‘" 


where 


e’ = e + e,d + e,d^ + e,d^ new constant term 

o o lo 2 o 3o 

ej = 6j + 202 d^ + 303d^ new coefficient of the 1st power of new 

time variable d ' 

Cj = 63 + 303d^ 


should be reduced by modulus 2 tt (= 6.283 185 307 179 586) using double pre- 
cision arithmetic when angular units of radians are used. For short periods of 
time , the ccntributions of the d ' and d ' terms are usually insignificant, so 
that the expression for the calculation of any mean element reduces to the sim- 
ple form 


e 

m 


e;d- 


d' = days elapsed from 
new epoch. 
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In most cases, single precision arithmetic may be used for a period of one 
year from the new epoch without significant loss of accuracy. The fastest mov- 
ing body for which mean elements will normally be required for current space 
work is the Moon, whose mean orbital motion is approximately 13.1 degrees per 
day. Letting e'^ be the mean motion of the Moon, the error associated with the use 
of single precision arithmetic may be calculated as follows: 

(e| ± Ae|) X (d' + Ad') = e| d' ± d' Aej ± e| Ad' 


neglecting the Ae| Ad' term. 

For single precision accuracy (approximately 8 decimal digits on the IBM 7094), 
Ae^ = ± (1 X 10“®) degrees/day for the mean motion of the Moon 


Ad' = ± (1 X 10 ^) days for d' = 1-9 days 

= ± (1 X 10“®) 10- 99 " 

= ± (1 X 10“®) 100-999 " 


SO that 


(e| ± Ae|) X (d' ± Ad') = e| d' + (1 x 10“®)° d’ ± 13?1 Ad' 


r 


(1 X 10“®) (9)'^ ± 13?1 (lx 10“ 0 




1 - 9 day 


= e|d' + ^(1 X 10“®) (99)° ± 13?1 (lx 10“®) V 10- 99 " 


^1 X 10“®) (999)° ± 13?1 (lx 10“®) 


100-999 
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10.30 X 10"6 


r 


9f00 X 10"® ± 1?3 X 10~® 


= e;d' ± < 9?90 X 10'5 + 1?3 X 10"5 = 


9?99 X 10"'’ ± 1?3 X 10'^ 

V- 


e[d' ± ^ 11.20 X 10“S 

J1.22 X 10"» 


/* 


= e;d'±^ ':s 


■’05 

for 1- 9 days 

■JS 

" 10- 99 " 

S'.'O 

" 100-999 " 


Tables of the mean elements of the Sun and Moon appear at the end of this 
report. These mean elements are listed in 6 combinations of units: 

Angular Units : Degrees, radians, revolutions 
Time Units : Days, Julian Centuries 


The tables are derived strictly from the expressions listed in the American 
Ephemeris and The Explanatory Supplement. These two references in turn list 
the elements as given by Newcomb^’'^ for the theory of the Sun and by Brown^®' 


for the theory of the Moon. In each case, the mean elements listed by all of 
these volumes are in units which are awkward and inefficient for computer 
usage. 


The enclosed tables use the same epoch, January 0.5, 1900, as used by the 
American Ephemeris and The Explanatory Supplement. Sufficient accuracy is 
carried in the coefficients listed in the tables to maintain a numerical accuracy 
of "01 over a time span of one century from the epoch. 


The above method of updating mean elements was programmed for the IBM 
7094 computer. In order to compare the method with the formulae used ly the 
Jet Propulsion Laboratory in their trajectory programs, new coefficients were 
computed for the epoch January 1.0, 1950. Significant discrepancies were found 

in the following 3 elements:^’h(io) 


(^)Newcomb, “Astronomical Papers Prepared for the use of the American Ephemeris,” Vol. 6, Part 1. 

Brown, E.W., “Tables of the Motion of the Moon." 
f^^Kalensher, B.E., "Selenographic Coordinates,” JPL Technical Report No. 32-41, page 23. 
(lO)Holdridge, D.B., “Space Trajectories Program for the IBM 7090,” JPL Technical Report 32-223, 
page 69. 


11 



Moon — Mean Anomaly 


(1) 

215?531 46 + 

13?064 

993 

d 

(Sandifer) 

(2) 

21 5?540 13 + 

13?064 

992 

d 

(JPL) 


+ f008 67 - 

?000 

001 

d 

difference (2-1) 

Moon - Argument of Perige6(’’1 





(1) 

196?731 199 + 

0?i64 

357 7 

d 

(Sandifer) 

(2) 

196?745 632 + 

0?164 

358 6 

d 

(JPL) 


+ ?014 433 + 

?000 

000 9 

d 

difference (2-1) 

Sun — Mean Anomaly d ^ ) 





(la) 

358?000 682^’^’ + 

0.985 

600 3 

d 

(Sandifer) 

(lb) 

358.000 670^’^' + 

0.985 

600 3 

d 

(Sandifer) 

(2) 

358?009 067 + 

0.985 

600 5 

d 

(JPL) 


+ ?008 397 + 

?000 

000 2 

d 

difference (2- la) 


Minor discrepancies were found in the other elements listed by JPL. The ele- 
ments listed above are used by JPL in the calculation of the lunar libration terms 
for conversion to Selenographic coordinates. The discrepancies noted are of the 
order of 1 x 10~^ degrees in the constant term or 2 x lO"** radians. A further 
source of discrepancy will be encountered when using the JPL formulae for the 
year 1966. During this year, the value of d is approximately 365 days x 16 years 
or 5 X 10^ days. Since the coefficients of the d terms carried by JPL in the 
above 3 cases is of the order of 10"® to 10"^ degrees, the error accumulated 
after 5 x lO^ days is 5 x lO"^ degrees or approximately 10"'* to 10~^ radians. 

It is therefore possible that a transformation to a Selenographic coordinate sys- 
tem using the JPL values for the mean elements could reduce the number of 


is also on internal inconsistency in JPL Technicol Report No. 32-223 on poge 69. The 
mean anomalies of the moon ond sun are defined as g = « - F* and g * = L - P and the argu- 
ment of perigee of the moon is defined as O) = F’ — Q. However, if the numerical values of 
«, L, F, F' and Q listed on page 68 are substituted into these equations, the results are not 
the same as the numerical expressions for g, g ' and o) listed on page 69. 

Corresponds to the value 358° 28' 33'.'04 updated from epoch Jan. 0.5, 1900. 

(13) " " " " 358° 28' 33'.'00 updated from epoch Jan. 0.5, 1900. 
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significant figures in a position vector from 8 significant digits to 3-4 significant 
digits, due to these causes only. As a double check, both the JPL formulae and 
the formulae computed ty the 7094 pr(^ram were used to calculate the 3 elements 
for a date in 1966, using double precision arithmetic in both cases. These values 
were then compared with the values tabulated in the American Ephemeris. In 
each case the JPL values were off by the same order of magnitude noted above, 
whereas the values predicted by the method described in this report, using a 
January 1.0, 1950 epoch, agreed exactly with the tabulated values.* It is therefore 
recommended that the expressions for the mean elements for the Sun and Moon, 
when required, be taken from the tables of this report. 

These above mentioned discrepancies are significant in that they illustrate 
a very bad computing procedure commonly used by many reputable organizations 
and individuals. There exist today a proliferation of epochs which have been chosen 
for various reasons. One of the reasons commonly advanced is that the use of dou- 
ble precision arithmetic is avoided by using an epoch close to the computing date. 
However, the method normally employed is to precompute the updated coefficients 
of the mean elements of the Sun and Moon on a desk calculator, and then enter 
these updated coefficients as constants in the computer. This method invariably 
leads to the kind of confusion encoimtered in the referenced JPL reports. This 
error is usually compounded in that other organizations and individuals adopt 
these erroneous values as published and thus perpetuate the error. 

As a case in illustration, the IBM group currently engaged in preparing the 
AIMP real-time program are using values of the coefficients of the mean ele- 
ments of the Sun and Moon which appeared in a report by Carson and Clark of the 
Data Operations Branch. These values, in turn, appear to have been adopted by 
Carson and Clark in good faith from the previously referenced JPL reports. 

In addition, the argument advanced — that the use of double precision arithmetic 
is avoided with the use of precomputed updated coefficients — is quickly outdated 
in a few years when the coefficients of the first power of time multiply into an 
increasingly larger nmnber. 

It would be preferable to have the various national observatories responsi- 
ble for maintaining ephemerides publish the elements in a form suitable for 
computers. If it is desired to minimize the use of double precision arithmetic, 
it is a fairly routine and simple procedure to program the computer to update the 
coefficients automatically and internally, thereby avoiding the type of errors en- 
countered in precomputing the coefficients on a desk calculator. 


13 



CALCULATION OF THE SELENOCENTRIC POSITION VECTORS OF 
THE EARTH AND SUN NEGLECTING PERIODIC PERTURBATIONS 


Initialization 


1 . Choose the Julian Date of 4th stage injection ( = JD_!, ) as the new epoch to 

replace the epoch JD^' = 2415020. 0 for the calculation of the mean orbital 
elements of the Earth and Sun.* (See note below.) 


2. Compute 


A. djj = JDq - JDjj in double precision. (JD^ must be stored as a double 
precision word. JD^^ may be stored as a single precision word. ) 

B. d^^, d^^ in single precision. 


3. Compute for JD^ the numerical values of the following mean elements: 


^3 (Tables le, 3) 

II 

S3 


Reduce all angles by modulus 2v using double precision arithmetic {2v - 
6.283 185 307 179 586) 

4. Compute the matrices 

R' = R 3 (“^ 3 )Ri(“i' 3 ) (Table 3 for i’3 , sin i'j , cos i'j ) 

R" = R 3 (u 3 )Rj(i'^) R 3 (Q’ 3 ) (NOTATION for definitions of Rj, R 3 ) 

5. Compute the coefficients of the sin ng and cos ng terms of the expansions 
of C^ and , using the algorithm in Table 4, to the 4th power of . 


*Any convenient epoch may be chosen. Time of 4th stage injection is chosen here os one pos- 
sible convenient epoch for AIMP. 
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6. 


Compute the new numerical coefficients of the expressions for , g ’3 , , 

and valid for the new epoch JD^ , according to the method described in 
the previous section. 

Calculation For Any Time JD After Initialization 


1 . Compute 

d ' = JD - JD^ = number of days elapsed since 4th stage injection (or 
other new epoch) 

2 . Compute Geocentric position vector of the Sun, referred to the ecliptic and 
mean equinox of date 


Compute: 

a. g^; reduce by modulus 2n 

b. C^, using algorithm of Table 4 

c. u ' = o)' + g' + C' ; reduce by modulus 2tt 

o o o ’ 

cos U * 

o 

sin Geocentric position vector of the Sun, ecliptic and 

mean equinox of date. 

_ 0 _ 

3. Compute Selenocentric position vector of the Earth, referred to the ecliptic 
and mean equinox of date. 



Compute: 

a. g' 3 ; reduce by modulus 2 tt 

b. C'^ , r’^ (Table 3)* 

c. U 3 = cJ'^ + g '3 + C* 3 *; reduce by modulus 2tt 
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d. 




COS u 


sin u. 


II 


r'j = R'q'3 Selenocentric position vector of the Earth, referred to 
the ecliptic and mean equinox of date 


4 . Compute Selenocentric position vector of the Sun, referred to the ecliptic 
and mean equinox of date. 


Sun 

7" =7" +7' Moon 

o 3 o 

5. Transform 7 '^ and 73 from ecliptic 

nate system formed by using the mean Earth-Moon orbit plane as the funda- 
mental plane and the mean Moon to Earth line at JD^ as the X axis (JD^' = 
4 th stage injection). 

7." ^ 

o o 

7 ; - R"7’’ 



and mean equinox of date to a coordi- 



TABLE la 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 


CENTER 
X-Y PLANE 
X AXIS 


Units 


ANGULAR 

TIME 


Numerical Values 


Geocentric 
Ecliptic of date 
Mean equinox of date 


* = revolutions, ° = degrees, ' = minutes, " = seconds 
T = number of Julian centuries elapsed since 12 hours ephemeris time, 
January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 


L' 


r' 

o 

n* 


- 

279°41'48':04 + 


279°41‘27"54 + 

= + 180° 3 

99°41'27':54 + 


281°13’15';00 + 


0° 

1 

f358°28' 33'.‘00^ 

- L‘ - r’ - ^ 

1 

O O j 

|358°28' 33104 

- r‘ - O' 

o o 


- r’ 

281° 13’ 15': 00 + 

o 


= l ; - 


= K 

279°41‘48':04 + 



1°43’ 09'.'03 T + 1'.’63 + '.'012 


1’.’541 


+ 99*'359"02’59:’10 T - 


- '.'012 


23°27' 08'.'26 - 
.016 751 04 - 


^'540 

1°43‘09'.'03 T + 1'.'63 + '.'012 

[^1.089 
|l'.'090 T* 

46'.'845T - '.'0059 + '.'00181 

.000 041 80 T - .000 000 126 T* 


‘This is eouivolent to T„' = 18^ 38*^ 45?836 + lOO'^ 00^ 03^ 0^542 T 4 ?0929 

• w 

^Constant term according to Newcomb is 358 °28' 33.0. 
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TABLE la (Continued) 


L'„, t’ 100'* 0°46'08'.'13 

^ o 

r' l°43'09r03 

o 

g' 99’*359°02'59'.'10 


129 602 768'.' 13 = 8 640 184?542 

6 189'.' 03 
129 596 579'.' 10 


Conversion of r'^to ', 


Constant Term 


T Term 


18" 

38" 

45 ® 

?836 

18" 38" 45?836 


270" 

9° 30' 

11' 15" 
12'.' 54 

279° 41' 27'.' 54 


100 " 

00 " 


03" 

04?542 

100" 00" 03" 04?,542 


100 ** 

00 ° 

0° 45' 

1' 08 '.' 13 

100" 00° 46'08'.'13 
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TABLE lb 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 


CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 

ANGULAR : ° = degrees 

TIME ; T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 

Numerical Values 


L' 


279?696 677 778 + 36 000?768 925 000 T 


?000 302 500 
?000 302 778 


T ' 
o 


r: 

279?690 

983 

333 

+ 

36 

000?768 

925 000 

T 

+ 

?000 

387 

083 

T2 



(r 

__ o 

- t’ + 180° 

O 

z: 

99?690 

983 

333 

+ 

36 

000?768 

925 000 

T 

+ 

?000 

387 

083 

T2 



r* 

o 



281?220 

833 

333 

+ 


1?719 

175 000 

T 

+ 

?000 

452 

778 

T2 

+ 

?000 003 333 

Q’ 

o 



0?0 



























r 









358?475 

833 

333 







?000 

150 

278 

j2 



g'o 

= L* - r' 

o o 

H 

^358?475 

844 

444 

+ 

35 

999?049 

750 000 

T 

1 

JOOO 

150 

000 

j2 


?000 003 333 

(l) * 

o 

^ r* - O' 

o o 


















^ r* 

o 


281?220 

833 

333 

+ 


1?719 

175 000 

T 

+ 

?000 

452 

778 

j2 

+ 

fOOO 003 333 

K 

- l; - 










1 

f?000 

302 

500 

T2 






279?696 

677 

778 

+ 

36 

000?768 

925 000 

T 

+ ^ 

1 

j^?000 

302 

778 

^2 



e 

m 



23?452 

294 

444 

- 


?013 

012 500 

T 

- 

?000 

001 

639 

j2 

4 

?000 000 503 

e' 

o 



.016 

751 

04 

- 


.000 

041 80 

T 

- 

.000 

000 

126 

T2 
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TABLE Ic 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate ^stem 


CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 


ANGULAR 

TIME 


” = degrees 

d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days exactly. 


Numerical Values 



1 

"?226 749 X 10'® d® 

l; 

o 

279?696 677 778 f ?985 647 335 387 d + ^ 



1 

?226 957 X 10-'® d® 

'^o 

279?690 983 333 i ?985 647 335 387 d + 

?290 151 X 10'® d® 

= r'+180" 

o o 

99?690 983 333 + ?985 647 335 387 d i 

?290 151 X 10'® d® 

r' 

o 

281?220 833 333 + ?000 047 068 446 d + 

?339 394 X 10'® d® 

.0’ 

o 

0?0 


1 

f358?475 833 333 

f?112 646 X 10' d^ 

K = ‘'o-' o = •{ 

+ ?985 600 266 940 d - J 

1 

1 

|358?475 844 444 ] 

1 ?112 437 X 10' '® d® 

of - r'4i‘ 



o o o 



: r” 

281?220 833 333 + ?000 047 068 446 d + 

?339 394 X lo'® d® 


?068 408 X 10'® d® 


o o r> 

L* 


279?696 677 778 + ?985 647 335 387 d + 


23?452 294 444 - ?000 000 356 263 


?226 749 X 10'*^ d® 

?226 957 X 10' d® 
d - ?001 228 X lO"'® d^ + ?010 318 x 10''® d® 


.016 751 04 - .000 000 001 144 421 d - .000 094 447 x 10' d® 
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TABLE Id 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 


CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 

ANGULAR ; '' = radians 

TIME : T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 


L’ 


= 4^881 627 934 112 + 628T331 950 990 910 T + 


Numerical Values 

"sr 279 621 X 10 ' « 

Sr 284 469 X 10' 6 T’ 

t' = 4 f 881 528 547 307 + 628 f 331 950 990 910 T + 6 f 755 879 x 10 " ® 

o 

= ’■'o + ^' = 1-739 935 893 717 628 f 331 950 990 910 T 6 f 755 879 x 10‘6 T^J 

r 030 005 264 168 T 7 T 902 463 x 10 ‘ ® + T 058 178 x 10 ' » 


r” 

o 

n* 

o 

g; - K - r; 

r-i* r r' - n' 

o o o 

= r' 


= 4f908 229 466 869 

= oro 

|6r256 583 580 497 
16:256 583 774 423 


628:301 945 726 742 T - 


2:622 842 x 10' •> 
2f617 994 X 10® 


- :058 178 X 10' « 


4:908 229 466 869 + :030 005 264 168 T + 7:902 463 x 10' « :058 178 x 10' « 


L* - O' 










r 



j2 

O O 










5^279 

621 

X 10'« 

K 

= 4r881 

627 

934 

112 

+ 628^331 

950 

990 

910 

T + J 

1 














1 5^284 

469 

X 10' « 

j2 















= f409 

319 

755 

203 

rooo 

227 

110 

969 

T - 

r028 

604 

X 10'® 

T2 


= .016 

751 

04 


.000 

041 

80 


T - 

.000 

000 

126 

T2 


:008 775 X 10' « 
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TABLE le 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 


CENTER 
X-Y PLANE 
X AXIS 

Units 

ANGULAR 

TIME 

Numerical Values 


Geocentric 
Ecliptic of date 
Mean equinox of date 


: ' = radians 

: d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days exactly. 


r003 957 513 x 10' d^ 

r003 961 147 X 10 ' ‘2 

non v in' * 2 2 


L* 

0 


- 4:881 

627 

934 

112 

f 

f017 

202 

791 

266 

007 

7 ' 

0 


- 4f881 

528 

547 

307 

+ 

f0l7 

202 

791 

266 

007 


1 ' 

0 

71^ If 739 

935 

893 

717 


f017 

202 

791 

266 

007 

I” 

0 


4f908 

229 

466 

869 

+ 

fOOO 

000 

821 

499 

361 

O' 

0 


- OfO 












f6f256 

583 

580 

497 







Ky 

l;^ 

< 




+ 

f017 

201 

969 766 

646 



|6f256 

583 

774 

423 







0 

f 

0 

- O' 

0 











i*' 

0 

4f908 

229 

466 

869 

+ 

fOOO 

000 

821 

499 

361 


K 

- O' 

0 











L,; 

4f881 

627 

934 

112 

\ 

r017 

202 

791 

266 

007 

m 


f409 

319 

755 

203 

- 

fOOO 

000 

006 

217 

959 

0 


.016 

751 

04 


- 

.000 

000 

001 

144 

421 


rOOl 966 037 X 10" 


- iri93 948 X 10‘2‘ d’ 


I roOl 962 403 X 10 ' ' 2 d^ 


|r003 957 513 x lo'^ d^ 
lr0O3 961 147 X 10' *2 d^ 


12 ,j 2 
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TABLE If 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 

CENTER : Geocentric 

X-Y PLANE : Ecliptic of date 

X AXIS : Mean equinox of date 

Units 

ANGULAR : * = revolutions 

TIME : T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 




I ?840 278 X 10‘« 

?776 935 216 049 + 100*002 135 902 778 T + 

|^?841 049 X 10 s T* 

?065 145 262 346 - ?000 036 145 833 T - ?004 552 x 10' « 4 1?397 x 10 ® 

.016 751 04 - .000 041 80 T - .000 000 126 
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TABLE Ig 

MEAN ORBITAL ELEMENTS OF THE SUN 


Coordinate System 


CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 


ANGULAR 

TIME 


* = revolutions 

d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days exactly. 


Numerical Values 


L' 


O' = i' t 







f?629 

858 X 

10- >5 

?002 

737 

909 264 

963 d 

A 








|?630 

436 X 

10- ‘5 

?002 

737 

909 

264 

963 d 

+ ?805 

975 X 

10- '5 

CM 

o 

o 

oc . 

737 

909 

264 

963 d 

+ ?805 

975 X 

10-' = 

?000 

000 

130 

745 

684 d 

+ ?942 

762 X 

10"* = 


L' 


= r - n' 


?o 

766 203 704 P'312 905 x 10' 

+ ?002 737 778 519 279 d _< - ?190 023 x 10'^* d^ 

?995 766 234 568 ?312 326 x lO''’ 


?781 168 981 481 + ?000 000 130 745 684 d + ?942 762 x 10' d^ i ?190 023 x 10'^' d^ 


F' . 1.' - O' 

o o o 

L' 











1 

r?629 

858 X 

10' 

1 5 

d2 

?776 

935 

216 

049 + 

?002 

737 

909 

264 

963 

d +J 

1 














1 

1^630 

436 X 

10' 

1 5 

d2 

?065 

145 

262 

346 - 

?000 

000 

000 

989 

619 

d - 

?003 

412 X 

10' 

1 5 

d2 

.016 

751 

04 


.000 

000 

001 

144 

421 

d - 

.094 

447 X 

10' 

1 5 

d^ 
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TABLE 2a 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER Geocentric 

X-Y PLANE : Ecliptic of date 

X AXIS : Mean equinox of date 

Units 

ANGULAR : ^ = revolutions, "" = degrees, ’ = minutes, " = seconds 

TIME : T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 

Numerical Values 


L' 

r 

270"26‘02'.'99 

-f 

1336*307°52'59'.'3l T 

- 

4" 08 

+ 

:0068 

!’• 

r 

334"19 46:40 

+ 

11'*109°02'02'.'52 T 

- 

37 '.'17 

- 

'.'045 

a* 

= 

259°10’59:79 

- 

5'*134°08'31'.'23 T 

+ 

7 '.'48 

+ 

'.'008 

g' 

= V -I ' 

295°06*16:59 

+ 

1325'*198°50'56'.'79 T 

+ 

33'.' 09 

+ 

'.'0518 

Cl)* 

r* -0‘ 

75°08’46V61 

+ 

16'*243‘’10'33"75 T 

- 

44'.'65 

- 

'.'053 

CjJ* 

f 180" - 

255°08’46:61 


16*243°10’33'.’75 T 

- 

44’.' 65 

- 

'.'053 

F 

= L’ ~n* 

11 "15‘ 03:20 

+ 

1342'* 82°0r3O'.'S4 T 

- 

11'.' 56 

- 

'.’0012 

0 * 

V -O' + 180" = 

191"15' 03:20 

+ 

1342'* 82°0r30'.'54 T 

- 

11'; 56 

- 

'.'0012 

D' 

-- L' -L' 

350"44'14:95 

4 

1236'*307°06’51"18 T 

_ 

5'; 17 

+ 

';0068 T ’ 


i' = 5'’08‘43V428 

I' - 1°32'06" 

e' :r .054 900 489 

a’ = 60.2665 Earth radii 

n 384,399.3537 km (using 1 Earth radius = 6378.3255 km.) 


sin i' 


+.089 683 453 

cos i' 

= 

i . 995 970 322 


sin I' 

r 

4.026 

787 

60 

cos I' 


4.999 

641 

15 
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TABLE 2b 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER 
X-Y PLANE 
X AXIS 


Units 


ANGULAR 

TIME 


Numerical Values 


Geocentric 
Ecliptic of date 
Mean equinox of date 


° = degrees 

T = number of Julian centuries elapsed since 12 hours ephemeris time, 
January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 


L‘ 


= 270?434 163 889 

+ 

481 2677883 

141 667 T - 

7001 

133 

333 

T2 + 

?000 001 889 

r’ 


n 334?329 555 556 

+ 

4 0697034 

033 333 T - 

7010 

325 

000 

T2 _ 

?000 012 500 

n’ 


- 259?183 275 000 

- 

1 9347142 

008 333 T 4 

7002 

077 

778 

) 

‘.’000 002 222 

r' 

: L’ -!’• 

296V104 608 333 

+ 

477 1987849 

108 333 T + 

7009 

191 

667 

1 

?000 014 389 

(•/)' 

P' -Q' 

- 75? 146 280 556 

■1 

6 0037176 

041 667 T - 

7012 

402 

778 

T2 _ 

'.’000 014 722 

CD 

+ 180° 

- 2557146 280 556 

+ 

6 0037176 

041 667 T - 

7012 

402 

778 

T2 _ 

‘.’000 014 722 

F' 

- L' -0.' 

r: 117250 888 889 

+ 

483 2027025 

150 000 T - 

7003 

211 

111 

T2 _ 

'.’000 000 333 

ir 

= L' -n' +180° 

^ 1917250 888 889 

+ 

483 2027025 

150 000 T - 

7003 

211 

111 

T2 _ 

?000 000 333 

D' 

L' -L' 

o 

: 3507737 486 111 

-1 

445 2677114 

216 667 T - 

7001 

436 

111 

T2 f 

?000 001 889 

i* 


- 57145 396 667 










1?535 000 
.054 900 489 

60.2665 Earth radii 

384,399.3537 km. (using 1 Earth radius = 6378.3255 km.) 


sin i 
cos i ' 


f.089 683 453 
t.995 970 322 


sin r 

cos r 


(.026787 60 
1.999 641 15 
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TABLE 2c 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER 
X«Y PLANE 
X AXIS 

Units 

ANGULAR 

TIME 

Numerical Values 
V 

[■” 

Q’ 

g’ - V -r* 

w' + 180 ° Q 

F’ :r V -O' 

f'’ - V -fr 4 180 ° - 

D’ - L’ -L' 

o 

i' 

I’ 

e’ r 

a* = 

sin i' = 

cos i' ^ 

sin I* 

cos I* c 


Geocentric 
Ecliptic of date 
Mean equinox of date 


: ” = degrees 

: d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days exactly. 


270?434 163 889 + 13?176 396 526 806 d - ?849 527 x 10 " 12 ^2 + 

334T329 555 556 + ?H1 404 080 310 d - 7?739 441 x 10 " ^2 _ 

259?183 275 000 - ?052 953 922 199 d + 1?557 466 x 10" ^ 

296?104 608 333 4 13?064 992 446 498 d 4 6?889 914 x 10'^^ 4 

75?146 280 556 + ?164 358 002 510 d - 9?296 906 x 10' ^2 _ 

255?146 280 556 + ?164 358 002 510 d - 9?296 908 x 10 " ^2 _ 

11?250 888 889 + 13?229 350 449 008 d - 2?406 993 x lO'^^ ^2 _ 

191? 250 S88 889 4 13?229 350 449 OOR H - 2?406 993 x lo '^ - 

350?737 486 111 + 12?190 749 191 421 d- 1?076 484 x + 

5?145 396 667 
1?535 000 
.054 900 489 


?038 765 X 10“ d^ 
?256 531 X 10' ** d^ 
?045 605 X 10** d^ 
T295 295 x 10*» d^ 
?302 136 X 10' d^ 
?302 136 X 10' d^ 
?006 841 X 10** d^ 
?006 841 X 10'*® d® 
?038 765 X 10' '« d® 


60.2665 Earth radii 

384,399.3537 km. (using 1 Earth radius ~ 6378.3255 km.) 


+.089 683 453 
+.995 970 322 

+.026 787 60 
+.999641 15 
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TABLE 2d 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER : Geocentric 

X-Y PLANE : Ecliptic of date 

X AXIS : Mean equinox of date 

Units 

ANGULAR : " = radians 

TIME : T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 

Numerical Values 

L' ^ 4f719 966 569 739 4 8 399:709 144 925 395 T - rOOO 019 780 398 + f032 967 x 10* 

r" = 5:835 151 531 174 t 7U018 041 257 371 T - :000 180 205 245 - :218 166 x 10 * 

0 ' ^ 4:523 601 514 852 - 33:757 146 246 552 T t :000 036 264 063 + :038 785 x 10* 

e' - L' -T' ^ 5:i68 000 345 745 + 8 328:691 103 668 024 T + :000 160 424 847 i :251 133 x 10 * 

I " -O' - i:311 550 016 322 i 104:775 187 503 924 T - :000 216 469 309 - :256 951 x 10 * 

+7?" = 4:453 142 669 912 < 104:775 187 503 924 T - :000 216 469 309 - :256 951 x 10 * 

F' L' -O' = :i96 365 054 887 + 8 433:466 291 171 947 T - :000 056 044 462 - :005 818 x 10'* 

O' L' -O' +77^ " 3:337 957 708 477 + 8 433:466 291 171 947 T - :000 056 044 462 - :005 818 x 10'* 

D' L' -L' n 6:121 523 942 807 + 7 77U377 193 934 485 T - :000 025 064 867 + :032 967 x 10 * 

O 

i' ^ \089 804 114 

r = r026 790 804 

- .054 900 489 


a’ - 60.2665 Earth radii 

= 384,399.3537 km. (using 1 Earth radius “ 6378.3255 km.) 


sin 

i' 

- 

+ .089 683 4 53 


cos 

i' 

- 

t . 995 970 322 

- 

sin 

I' 


i. 026 787 60 


cos 

r 

□ 

+.999641 15 
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TABLE 2e 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate l^stem 

CENTER 
X-Y PLANE 
X AXIS 


Units 


angular 

TIME 


Numerical Values 
V 

r 

Q 

: L’ -r* 

: r -,a‘ 


Geocentric 
Ecliptic of date 
Mean equinox of date 


= radians 

d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0), One Julian century = 36525 days exactly. 


V -Q‘ 

L’ -n* +7/^ 

L’ -L’ 


4^719 966 569 739 + 
5^835 151 531 174 + 
4:523 601 514 852 - 
55168 000 345 745 + 
15311 550 016 322 + 
45453 142 669 912 + 
5196 365 054 887 + 
35337 957 708 477 + 
65121 523 942 807 + 
5089 804 114 
5026 790 804 

. 054 900 489 


60.2665 Earth radii 

384,399.3537 km. (using 1 Earth radius " 6378.3255 km.) 


941 

147 d - 

5014 

827 

044 X 

10-12 

001 

571 d - 

5135 

078 

733 X 

10-12 

294 

225 d 4 

5027 

182 

914 X 

10-12 

939 

576 d + 

5120 

251 

689 X 

10-12 

295 

795 d - 

5162 

261 

647 X 

10-12 

295 

795 d - 

5162 

261 

647 X 

10-12 

235 

372 d - 

5042 

009 

958 X 

10-12 


0*70 J 
0 / i VJ — 

*^042 

009 

958 ^ 

10-12 

675 

140 d - 

5018 

788 

191 X 

10-12 


3-2 1 ,J3 
3-21 ^3 
3-21 
3-2^ d^ 

3-2 1 ,j3 

3-21 

3-21 ,j3 
3-21 


Sin 1 
cos i' 


+.089 683 453 
+.995 970 322 


sin I' 

cos I' 


+ .026 787 60 
+ .999 641 15 
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TABLE 2f 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 

ANGULAR : ^ = revolutions, ° = degrees, * = minutes, ’ = seconds 

TIME : T = number of Julian centuries elapsed since 12 hours ephemeris time, 

January 0, 1900 (Julian Date = 2415020.0). One Julian century 
= 36525 days exactly. 

Numerical Values 


V 


r' 

Q' 

r' 


Cl) 


r<) 


F' 

f)' 

D' 

i' 

I’ 


L’ -r- 
r-n- 

L’ -n' 

L’ -fi’ 

L’ -1^0 = 


?751 206 010 802 + 
?928 693 209 877 + 
?719 953 541 667 - 
?822 512 800 926 + 
?208 739 668 210 + 
?708 739 668 210 + 
?031 252 469 136 + 
?531 252 469 136 + 
?974 270 794 753 + 
?014 292 769 
?004 263 889 
.054 900 489 


1 336?855 230 949 074 
11?302 872 314 815 
5?372 616 689 815 
1 325?552 358 634 259 
16?675 489 004 630 
16?675 489 004 630 
1 342?227 847 638 889 
1 342?227 847 638 889 
1 236?853 095 046 296 


T - 3?148 148 X 10* 
T - 28?680 556 x 10' « 
T + 5?771 605 X 10* 
T 4 25?532 407 x 10'« 
T - 34?452 161 x 10 * 
T - 34?452 161 x 10' * 
T - 8?919 753 x 10‘* 
T - 8?919 753 x 10* 
T - 3?989 198 x 10* 


+ 5?247 X 10’ 

- 34?722 X 10' ’ 

+ 6?173 X 10' ’ 

T2 t 39?969 X 10"’ 
T2 - 40?895 X 10"’ 

- 40?'895 X 10"’ 

- ?926 X 10"’ 

T* - ?926 X 10" ’ 

t 5?247 X 10"’ 


a 


60.2665 Earth radii 

384,399.3537 km. (using 1 Earth radius = 6378.3255 km.) 


sin i' 
cos i ' 


+.089 683 453 
f.995 970 322 


sin I' 

cos r 


^.026 787 60 

+.999 641 15 
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TABLE 2g 

MEAN ORBITAL ELEMENTS OF THE MOON 


Coordinate System 

CENTER 
X-Y PLANE 
X AXIS 


Geocentric 
Ecliptic of date 
Mean equinox of date 


Units 


ANGULAR : = revolutions 

TIME : d = number of days elapsed since 12 hours ephemeris time, January 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days 

Numerical Values 


L' 


L'-r 

r- -Q’ 

+?5 

L' -O' 

L' -O' +?5 
L' -L' 


?751 206 010 802 + ?036 601 101 463 356 d - 2?359 797 x 

1*928 693 209 877 + ?000 309 4 55 778 640 d - 21?498 448 x 

?719 953 541 667 - ?000 147 094 228 332 d + 4?326 295 x 

?822 512 800 926 + ?036 291 645 684 716 d + 19?138 651 x 

?208 739 668 210 + ?000 456 550 006 971 d - 25?824 743 x 

?708 739 668 210 + ?000 456 550 006 971 d - 25?824 743 x 

?031 252 469 136 + ?036 748 195 691 688 d - 6?686 093 x 

?531 252 469 136 + ?036 748 195 691 588 d - 6?6S6 093 x 

?974 270 794 753 + ?033 863 192 198 393 d - 2?990 233 x 

?014 292 769 
?004 263 889 
. 054 900 489 


10“ 

1 5 

d2 

10' 

15 

d2 

10" 

15 

d2 

10" 

1 5 

d2 

10" 

15 

d2 

10' 

1 5 

d2 

10' 

1 5 

d2 

10' 

I 5 

d2 

10' 

15 

d2 


?107 680 X 
?712 585 X 

?126 682 X 

?820 265 X 
?839 267 X 
?839 267 X 
?019 002 X 
?019 002 X 
?107 680 X 


a 


60.2665 Eai th radii 

384,399.3537 km. (using I Earth radius 


6378.3255 km.) 


sin i' 


+ .089 

683 453 

cos i’ 


+ .995 

970 322 


sin I' 

= 

+ .026 

787 

60 

cos I' 

- 

+ .999 

641 

15 


0 , 

exactly. 

10'^* d^ 
10'^' d^ 
10’^* d^ 
10"^* d^ 
10-21 
10-21 ^3 
10-21 d2 

10-21 tU 

10-21 fp 
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TABLE 3 

MEAN ORBITAL ELEMENTS OF THE EARTH* 


Coordinate System 

CENTER 
X-Y PLANE 
X AXIS 


Units 


ANGULAR 

TIME 


Numerical Values 


Selenocentric 
Ecliptic of date 
Mean equinox of date 


' = radians 

d = number of days elapsed since 12 hours ephemeris time, January 0, 

1900 (Julian Date = 2415020.0). One Julian century = 36525 days exactly. 



: Q’ 

- 4^523 

601 

514 

852 - 

rooo 

924 

220 

294 

225 

d + 

ro27 

182 

914 X 

10" 

12 

+ 

T795 

965 X 

10‘ 

21 

d3 

g3 = 

^ g‘ 

- 5^168 

000 

345 

745 + 

T228 

027 

134 

939 

576 

d + 

ri20 

251 

689 X 

10' 

12 

+ 

5fl53 

876 X 

10' 

2 1 

d^ 

'^"3 '■ 

- C<>' + 77 '’ 

- 4T453 

142 

669 

912 + 

roo 2 

868 

588 

295 

795 

d - 

ri62 

261 

647 X 

10' 

12 

- 

5^273 

271 X 

10' 

21 

d3 

F" : 


^ 3f337 

957 

708 

477 4 

T230 

895 

723 

235 

372 

d - 

r042 

009 

958 X 

10' 

12 

d2 - 

rii9 

395 X 

10' 

21 

d^ 


COS I , - cos 1 


r089 804 114 
.054 900 489 

60.2665 Earth radii 

384,399.3537 km. (using 1 Earth radius - 6378.3255) 
+ .089 683 453 
+ .995 970 322 




001 507 032 - 


054 

838 

449 

cos 

e'i 

C 3 = 

ri09 

759 

636 

sin 

R3 

001 

504 

006 

cos 

2^3 

f 

roo3 

763 

418 

sin 

2R3 

000 

061 

878 

cos 

^63 

+ 

rooo 

178 

928 

s in 

3^3 

000 

003 

017 

cos 


+ 

rooo 

009 

721 

sin 

'*63 

000 

000 

161 

cos 

Se's 


rooo 

000 

568 

sin 

5^3 

000 

000 

009 

cos 


+ 

rooo 

000 

035 

s in 

6R3 

000 

000 

001 

cos 


+ 

rooo 

000 

002 

s in 

7gj 






f'3 ' e'j t 

C’a 










U 3 = Fj + - 

^3 + 

^^3 





*This table derived from Table 2e. 
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TABLE 4 


EXPANSIONS IN ELLIPTIC MOTION OF THE EQUATION OF THE 
CENTER AND THE MAGNITUDE OF THE RADIUS VECTOR* 


C = + 


« 1 , 5 , 107 , 

2 e - e + e | sin g 


[ 4 

[ 1 

103 
96 

1097 


5 0 _ 11 e" + 

4 e -24 ® + 192 ® 


3 43 5 95 7 

»3 _ + p7 

64 ® 512 ® 


4 451 , 

p4 _ „6 

480 


, 5957 , 

~ e' 

^ A f r\r\ 


960 4608 

1223 , 



960 

47273 , 


32256 


] 

j sin 2g 
J sin 3g 
j sin 4g 
j sin 5g 
J sin 6g 
j sin 7g 


r - a 




e + -g e3 - 


4-e2+ ^e^- 


3 

3 , 45 

8 128 


5 

® ^ 9216 


16 




COS g 


5 567 , 

e' 


5120 


3 

125 

384 

27 

80 


1 4 2 . 

+ ~ e” 

3 


e® + 


4375 

9216 


16807 , 

• 


46080 


J COS 2g 
J cos 3g 
j cos 4g 
J cos 5g 
J cos 6g 
j cos 7g|* 


‘'For example, refer to “Methods of Celestial Mechanics,” page 76-77, by Brouwer and Clemence. 
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